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Abstract 

Using a holographic method, we further investigate the relaxation towards the hy¬ 
drodynamic regime of a boost-invariant non-Abelian plasma taken out-of-equilibrium. 
In the dual description, the system is driven out-of-equilibrium by boundary sourcing, 
a deformation of the boundary metric, as proposed by Chester and Yaffe. The effects 
of several deformation profiles on the bulk geometry are investigated by the analysis of 
the corresponding solutions of the Einstein equations. The time of restoration of the 
hydrodynamic regime is investigated: setting the effective temperature of the system at 
the end of the boundary quenching to (r*) = 500 MeV, the hydrodynamic regime is 
reached after a lapse of time of 0(1 fm/c). 
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1 Introduction 


The possibility of employing the gauge/gravity duality correspondence 01 to the anal¬ 
ysis of far-from-equilibrium processes in strongly coupled systems has enlarged the 
realm of nonperturbative phenomena to which this theoretical method can be applied 
in a quite controllable way, while traditional approaches are less effective. Among the 
different systems, of particular interest is the one produced in ultrarelativistic heavy ion 
collisions (HI), as those taking place at RHIC and at LHC. The features of this system, 
for time scales larger than about 1 fm/c, seem to be well reproduced in a hydrodynamic 
setup involving a strongly coupled/low viscosity fluid [^[^, a framework allowing to cor¬ 
rectly describe experimental observables such as the light hadron transverse momentum 
spectra. In the hydrodynamic scheme of crucial importance are, after the pre-equilibrium 
regime following the collisions, the conditions at the time when the hydrodynamic behav¬ 
ior sets up. In particular, at this time the system stress-energy tensor Tif is important. 
We write it as 


Ar2 

Tl^ = ^diag{-e,p±,p±,p\\) , (1) 

in terms of the system energy density e and of p±,p\\ the pressures along one of the two 
transverse directions (with respect to the collision axis) and in the longitudinal direction, 
respectively]^ For example, in approaches based on the idea of an initial state described 
by color glass condensate with a saturation scale Qs of the order of a few GeV, the 
initial stress-energy tensor comes from classical chromo-electric and chromo-magnetic 
fields and has the form Ty oc diag(—e,e,e,—e). However, such an initial condition has 
been shown to be unstable, as soon as one-loop corrections in the strong coupling constant 
are switched on [^. 

In the analysis presented in our study, a prominent role is therefore played by the 
stress-energy tensor, which is needed to understand the features of the equilibrium 
regime, the time needed to reach equilibrium, and the properties of the dynamics driv¬ 
ing the system towards equilibrium. Questions to be answered are, for a system driven 
out-of-the-equilibrium, whether the late-time state is described by hydrodynamics, and 
whether the elapsed time needed to reach equilibrium is related to the dynamics during 
the out-of-equilibrium state. These issues concerning the transient process from a far- 
from-equilibrium to a hydrodynamic regime can be faced in the holographic approach , 
using as a guideline the conjectured correspondence between a supersymmetric, confor¬ 
mal field theory defined in a four dimensional (4D) space-time and a gravity theory in 
a higher dimensional space-time, a five dimensional Anti de-Sitter space (AdSs) times 
a five dimensional sphere [^. The gauge theory is defined on the 4D boundary of 
AdSs, and the connection has been extended to the case of nonlinear fluid dynamics [^. 
Non-equilibrium can be studied by solving in the bulk the 5D Einstein equations for the 
metric subject to suitable time-dependent boundary conditions. Information about the 


late-time regime can be obtained computing various invariants. For example, in Ref. 10 


the square of the Riemann tensor (the Kretschmann scalar) has been studied as an 
expansion in inverse powers of the proper time up to second order, finding that, in 
the asymptotic r —)• oo regime, the scalar 3?^ is free of singularities in the perfect fluid 
case. The boundary theory stress-energy tensor results from the solutions of the Einstein 
equations in the bulk. 


^For a review see [^. 

^ Along the paper, we refer to energy density and pressures without considering the factor ■ 
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Analyses in the case of boost-invariant fluid 11,12 , for some choices of the profile 


distortion, have shown the formation of a horizon in the bulk, and have given access to the 
components of the boundary theory stress-energy tensor. The investigations have been 
extended to samples of initial states characterized by different values of the components 
of the stress-energy tensor, without reference to the mechanism producing each state 


13-18 


As formulated in 11,12, the initial off-equilibrium states can be thought as being 


produced by external sources distorting the boundary metric for short time intervals. 
We focus on this issue, with a study of external distortion profiles, step-shaped or short- 
duration sequences repeating themselves with different intensities, with the aim of de¬ 
scribing phenomena where a small number of collisions takes place before the system 
starts evolving to thermal equilibrium. This study in the boost-invariant framework is 
useful for more involved cases where a lower symmetry is assumed for the system, for 
example in shock-wave collisions 


19-22 


There are several reasons to study complex distortion shapes in the boundary metric. 
The bulk Einstein equations are nonlinear, and different distorsions produce unrelated 
responses. The resulting effective temperature and entropy density follow different pro¬ 
files as a consequence of the boundary sourcing. Moreover, it is interesting to investigate 
whether the onset of the hydrodynamic regime is related to the distortion curve. Fi¬ 
nally, these analyses have to face problems concerning the stability of the solution of GR 
equations with different boundary conditions. 

The plan of the present paper, in which a few of the above issues are examined, is the 
following. In Section we describe how a system is taken out-of-equilibrium through the 
distortion of the boundary metric, as done in 11,12. We discuss a suitable procedure 


for solving the resulting 5D Einstein equations, providing a few details of the algorithm 
and the issue of the stability of the solutions. In Section we describe the profiles of 
the boundary geometry distortion we have chosen to study, with the motivations. In 
Section we discuss our findings, in particular concerning the equilibration time. The 
conclusions are presented in the last Section. Details dealing with the calculation of the 
boundary stress-energy tensor Tp are collected in appendix while in appendix we 
discuss some aspects of the numerical algorithm. 


2 Distorting the boundary 

The main idea, as in [^[^, is to study the effects of a time-dependent deformation 
of the metric of the 4D boundary. As a consequence of the deformation, gravitational 
radiation is produced and propagates in the 5D bulk, and a black hole is formed together 
with its horizon. In this way, one can trade the study of the reach of equilibrium in the 
4D boundary theory for the analysis of the time evolution of the geometry in the 5D 
bulk. 

We set the 4D coordinates = (x°, x^, x^, x^), identifying x^ = xy with the axis in 
the direction of collisions and along which the plasma expands. Boost-invariance along 
that axis is imposed, together with translation invariance and 0(2) rotation invariance 
in the orthogonal plane x_l = {x^, x^}. In terms of the proper time r and rapidity 
y, given by x*^ = rcoshy, xy = rsinhy, the 4D Minkowski line element is given by 
= —dr^ -|- dx^ -|- r^dy^. 

A distortion of the boundary metric, leaving the spatial three-volume unchanged 
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and respecting the imposed symmetries, can be obtained through a function 7 (t) which 
encodes information about a deformation: 

ds2 = ^ ^ . (2) 

The space-time with the metric Q is considered as the boundary of the 5D space-time 
in which the gravity dual is defined. We adopt 5D Eddington-Finkelstein coordinates, 
with r the extra-dimension coordinate and the boundary reached for r —>■ oo. In general, 
the 5D bulk metric can be chosen as 


ds‘^ = 2drdT — Adr"^ -|- Y,‘^e^dx‘]_ -|- ‘^^dy^ , (3) 


with the functions A, T, and B only depending on r and r to respect the adopted symme¬ 
tries. The behavior of these functions against the deformation allows us to describe how 
the bulk metric changes as a consequence of the distortion of the boundary. With the 


metric Q, at fixed x±_ and rapidity y, infalling radial null geodesics correspond to con- 

d'V 

stant r, while outgoing radial null geodesics are obtained from — = -^—. Therefore, 

dr 2 


for a generic function ^(r, r) the derivatives 


( 4 ) 

i = dri + ]^Adri ( 5 ) 


represent directional derivatives along the infalling radial null geodesics and the outgoing 
radial null geodesics, respectively. 

The 5D metric Q is the solution of Einstein’s equations with negative cosmological 
constant. In terms of ^(r, r), S(r, r) and B{r,T) the equations can be rephrased as [T^: 


E(S)' -h 2E'S - 2S2 = 0 

( 6 ) 

T.{B)' + ^ (i:'B + B't^ = 0 

( 7 ) 

A" + W’13 12 „ -k4 = 0 

( 8 ) 

S -F ^ - A't^ = 0 

( 9 ) 

E" + ^B'‘^T. = 0 . 

2 

(10) 


The five equations ([^-(10) can be considered as three dynamical and two 
equations. The condition that the metric Q produces the 4D metric Q at the 
for r — )■ oo, constrains the large r behavior of A{r,T), T,{r,T) and B(r,T): 


constraint 

boundary, 


r 

S(r, r) 
r 

B{r,T) 


— r 


0 




-> -xlogr-bqM • 

r^oo o 


( 11 ) 

( 12 ) 

(13) 
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The invariance of Q under the diffeomorphism r —)> r + A(r) has been exploited to 
impose the large r condition for A. 

We switch on the distortion of the boundary metric at the initial time t = Ti, starting 
from the AdSs bulk metric: 


ds'^ 


2drdT + 


—+ dx\ + 



(14) 


As discussed in taking the limit r —>■ oo and r —)• 0 with the Eddington-Finkelstein 
coordinates, Eq. ( |14[ ), is ambiguous. For this reason, we set Ti > 0. The initial conditions 
for A, E and B are therefore: 


Mr, n) 


Mr,n) 


B{r,Ti) 


Aini (t) 
^ini ) 


= r 


= r 



Bi, 




(15) 

(16) 

(17) 


For r —)• oo, the Einstein’s equations can be solved starting from the relations (|11[)- 
(13). As suggested in 11,12 , for A(r, r), E(r, r) and B{r,T) the large-r expansions can 
be written as 


Aasyir,T) = ^ [an(r) + On(r) logr]r^ , (18) 

n=0 

T,asyir,T) = ^ [Sn(T) +cr„(r) log r]r^"” , (19) 

n=0 

Basyir, r) = ^ [6n('r) + /3n(r) logr] . (20) 

n=0 


The conditions (0-([T3]) fix the n = 0 coefficients in (|18[)-([20l), as well as those for n = 1 
in Aasy The other coefficients can be determined imposing that Eqs. (|^-([^ are satisfied 
by the functions ([T^-([ 20 l), but this leaves two coefficients undetermined, 04 (r) and 64 (r). 


From Eq. (10) a relation follows between 04 and 64 : 


1296r®a4(r) + 1728r^a4(r) + G(r) 

if7no_4 / n I o_ 


( 21 ) 


with the primes denoting derivatives with respect to r, and the function G(t) expressed 
in terms of ■ 

G{t) = —576 + 3618r7'(r) — (69037'(r)^ + 14747"(r)) 

+ ^46087'(r)^ + 51997'(r)7''(T) + 4687 ^^^(r)^ 

+ 9187^(r)‘^ — 64807 ^(r)^ 7 ''(r) + 2167 ^'(r)^ — 5407'(r)7^^^(T) + 3067^^^(r)^ 

+ M ^32407'(r)^7''(r) + 1627 "(T) 7 ^^^(r) — 4597'(r)7^^^(T)^ . ( 22 ) 
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In our procedure, a 4 (T) and 64 (r) are determined by Eq. (21) and matching the 


computed functions S and B with their asymptotic expressions (19) and (20) for large r. 


Using 04 (r) and 64 (r), important observables can be obtained, such as the components 
of the boundary stress-energy tensor T/f, as discussed in next sections. 

The details of our method for solving Eqs. ([^-(10) are as follows. In the set of func¬ 
tions S, i], B, B and A, we treat a given function and its derivative ([^ as independent. 
At r = Ti, S(r, Ti), B{r, Ti) and A(r, Ti) are known, Eqs. To coimlete the set of 

initial conditions, we compute Ti{r,Ti) and B{r,Ti) solving Eqs. ^ and ([7), respectively. 

At r > Ti the algorithm is organized in steps. 

Using the definition ([^, we determine dr^{r,Ti) and drB{r,Ti), which allow us 
to obtain S and B at the subsequent time Ti + dr, if dr is small enough for the 
equations 

S(r, Ti -h dr) = I](r, n) -h dTdrT,{r, n) 


1 . 


B{r, Ti + dT) = B{r, Ti) + dTdrB{r, Ti) 


to be valid. 


2 . 


Once S(r, Tj -|- dT) and B{r,Ti + dT) are determined, we compute ai{Ti + dT) and 
64 (tj -t- dr). Indeed, 04 (r* -|- dT) can be obtained fitting the computed functions 


T,{r,Ti + dT) and B{r,Ti + dT) with their asymptotic expressions (19) and (20) in 
the range Vas ^ Tmax (later on in this section, it will be shown how to set Vas)- 
Then, we get 64 (rj -|- dT) from Eq. (21). 


3. 


i;(r, Ti + dT), B{r, Ti dT) and finally A{r, Ti + dT) are computed solving Eqs. 
(®; the boundary conditions are fixed using 
64 (t -h dr). 


( 20 ), (18), and 04 (ri -|- dr) and 


With the functions known at Ti + dT, the cycle is repeated until a chosen final value of r 
is reached. Typically, with the 7 functions considered in our analysis we use dr = 10“^. 

At each r step, our solution algorithm requires three values of the bulk coordinate 
'f'- fmax and Tmin, which are the maximum and minimum value of r used in the nu¬ 
merical calculation, and Vas, a value employed to determine 04 (r) from the asymptotic 
expressions (18)-(20) used in the range [ras,rmax]- Stability against variation of such 


parameters is used as a criterium to check the numerical results. Particularly relevant is 
the minimum value required to perform an excision in the bulk coordinates. This 


value is set in the following way. At fixed r, for each one of Eqs. (§-([T§ we construct 
the function TZi{r, r) defined as the absolute value of the l.h.s. of the equation divided by 
the sum of the absolute value of its addendi (see appendixfor definitions). At r = r*, 
we choose a small initial value of Tmin] then, at each r, is increased by steps of size 
typically 0.005, until the differential equations ([^-([^ and one constraint equation (10) 
are fulfilled, namely by requiring the condition TZi{rmimT) < 0.03 for i = 6,7,8,10. As 
a result, the bulk excision is always beyond the apparent horizon. In r-min, the function 
A{rmin,T) becomes negative but not large, hence possible singularities beyond the hori¬ 
zon are avoided. A similar procedure is carried out for We initially set r^s, and then 
we gradually increase it until the condition TZi (r^^, r) < 0.01 for i = 6,7, 8 ,10 is satisfied. 
Therefore, at this stage, we monitor the accuracy of the three differential equations and 
one constraint equation in and Vas- At the end of the procedure, we check that the 
condition TZi < 0.01 is fulfilled for all Eqs. ([^-( 10), i.e. f = 6 — 10, and in the whole (r, r) 
range. We find that this is indeed the case except for few tiny regions. In the largest 
part of the (r, r) plane deviations from zero are smaller than 0 ( 10 “^), as discussed in 
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appendix [B| During the numerical evaluation, no corrections are required to account for 
the finite value of the time step. In spite of the intricacy of the boundary conditions, 
stable solutions are found 0 


3 Models for the distortion of the boundary ge¬ 
ometry 


We are interested in investigating deformations of the boundary metric characterized by 
different duration, intensity and structure. We set the generic form for 7 (t): 


7 (r) = w 


Tank 


T -To 


N 

E' 

i=i 


+ > \lj{r,To,j) 


with 


and 


7i(D7-o,i) 


) = cs/j(T, ( 1 - 


= 1 - 


('T - 




(23) 


(24) 


(25) 


This expression generalizes the one used in , and represents sequences of N ” pulses”, 
each one of intensity cj and duration proportional to Aj, with the possibility of super¬ 
imposing a smooth deformation, of intensity w, which asymptotically changes the scales 
of the coordinates. After their generation, pulses travel in the longitudinal direction xy, 
driving the boundary state out-of-equilibrium. The response of the 5D bulk geometry 
to this deformation describes how equilibrium is reached at late times. We choose three 


different models. A, B and C, setting the various parameters in Eqs. (23)-(25). 


Model A represents two pulses in the boundary metric. The parameters are set to 
w = 0 and N = 2, ci = 1, Ai ^2 = I, t'o,! = f Ai, C 2 = 2. Moreover, two different values 
of To, 2 : To ,2 = and ro ,2 = |A 2 are considered, changing the time interval between 

the pulses. The final time of the distortion is = 5.25 and = 3.25 for the two cases, 
respectively, two quantities important for our discussion. 

In model B we set w = ^, rj = 1.2, tq = 0.25, A = 1, ci = 1, Ai = 1, ro,i = 4Ai, 
and the short pulse ends at r® = 5 while a tale of the distortion continues with r and 
approaches a constant value. 

Model C combines both the previous ones, and is obtained using 


w = -, ri = 1.2, To = 0.25 , 

5 

co,i = 0.5, Co ,2 = 1, Co ,3 = 0.35, co ,4 = 0.3, 

To,i = To-M.5, To, 2 = To-h 3.8, To, 3 = To6.1, ro,4 = To8.4, 
A, =0.8 (j = l,2,3,4). 


^ Other methods for the numerical solution of analogous General Relativity equations in the case of initial 


value problems are described in 23 and in references therein. 
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For this model the last short pulse ends at = 9.45. As mentioned, to avoid in the 
Eddington-Finkelstein coordinates the ambiguity in the limits r —)• oo and r —)• 0, in all 
cases the initial time is Tj > 0. The three profiles 7 (r) are depicted in Figs. [^[^and[^ 
The distortion profiles are chosen with the purpose of studying different situations. 
Model A is the repetition, in the boundary metric, of two short pulses of different inten¬ 
sity, Fig. In particular, we consider two cases with different interval of time between 
the pulses, ranging from the condition of distant perturbations to the case of overlap¬ 
ping pulses. This model can provide us with information on the horizon formation, in 
comparison with the case of isolated pulse studied in [^, and on the effects related 
to the time elapsed between the distortions. Model B has the purpose of studying a 
combination of two effects with very different time scales, a short pulse and a slow con¬ 
tinuum asymptotically producing a rescaling of the boundary coordinates. Model C, a 
combination of the previous cases, accounts for the effects of several pulses with var¬ 
ious intensities, and is closer to physical processes driving systems out-of equilibrium 
in relativistic heavy ion collisions. In all cases, the horizon formation and behavior are 
investigated, together with the various components the boundary stress-energy tensor, 
looking at the time when the hydrodynamic regime sets in. 


4 Results and discussions 


In the numerical calculation for the models previously described, the solutions of the 
bulk Einstein equations allow to obtain the apparent horizon and the event horizon, 
which define the system temperature T^ff. The apparent horizon is determined as the 
locus of points where the condition T,{rh,T) = 0 is fulfilled. The event horizon, which 
separates causally disconnected regions, is found drawing outgoing radial null geodesic 
curves and looking for the critical curve which separates the geodesics escaping towards 
the boundary from the ones plunging back in the bulk. Erom these curves, the function 
r/i(r) (r/i(r) being the position of the horizon at the proper time r) is determined, which 
defines the temperature Tejj(r). Asymptotically in proper time, the apparent horizon 
and the event horizon coincide. Erom the solution, S(r/i,r)^, the area of each horizon 
per unit of rapidity, can be computed. 

The three models share common features, which we discuss together with the differ¬ 
ences. Before presenting this analysis, let us recall a few results concerning the boundary 
stress-energy tensor Ty . 

As discussed in the Bjorken’s seminal paper [^, under the assumptions of homo¬ 
geneity, boost-invariance and invariance under rotations in the transverse plane with 
respect to the fluid velocity, the fluid stress-energy tensor has well-defined properties, 
and its components only depend on the proper time r, a manifestly boost-invariant con¬ 
dition. Imposing conservation and traceless conditions, the tensor components can 


^The apparent horizon is the outermost trapped null surface. To compute it, a foliation of spacetime at 


fixed times is considered. In each spacelike hypersurface, the apparent horizon is a closed surface such that 24 


H = 


Sb-K + S^S^Kab = 0 


with 7 ob the induced metric and Kab the extrinsic curvature of the hypersurface, K = g°‘^Kab, and s“ the 
outward-pointing spacelike unit vector, normal to the apparent horizon and tangent to the hypersurface. For 
the spacetime described by Eq. ^ this corresponds to the condition S = 0. 

^Studies of the behavior of the horizons in the gravity dual of a boost-invariant flow can be found in 25p6 





be expressed in terms of a single function /(r), that can be chosen to be —Tq : 

r^ = dza5r^-/(r), /(r) + ^r/'(r), /(r) + ^r/'(r), -/(T)-r/'(r)^ . (26) 

For a perfect fluid, the equation of state e = 3p, with p = p\\ = p±, fixes the r dependence: 
COTlst 

^('^) ~ stemming from the relations 

^4 /6 


P\\iT) = -e(r) -re'(r) 

(27) 

P±{t) = e(r) -P ^e'(T) . 

(28) 


Deviations from the ideal behavior can be included as corrections taking viscous effects 
into account, in a gradient expansion. On the gauge theory side, the gradient expansion 
corresponds to a late time expansion, with subleading terms identified with the dissipa¬ 
tive ones in the hydrodynamical theory. As shown in [10] , the function /(r) deviates 
from /(t) = const if the assumption of perfect fluid behavior is removed. For ex- 

1 . / N 1 , (s-1) 1 (2-s) 

ample, the case e(rj ~ — corresponds to py ^ "— * 


and p± 

" T“ 2t^ 

corrections to the perfect fluid behavior relate the energy density 

3.4. 


■■0 


Subleading 


e(r) = -TT Teffir) 


(29) 


to the definition of an effective temperature Tgff^r). The calculation in AA = 4 SYM 
gives 01 00: 


Teffir) = 


A 


(Ar):^/^ 


1 - 




-1 log 2 


67r(Ar)2/3 


-2127r2-p 51 log 2 - 24(log 2)2 
19447r3(AT)2 


This leads to 


e(r) = 


P\\{'r) = 


P±{t) = 


Svr^A'^ 


1 - 


2ci 


-P 


C2 


4(Ar)^/^ (Ar)2/3 (Ar)"^/^ 


vr^A^ 


1 - 


6ci 


-P 


5C2 


4(Ar)^/^ (Ar)2/3 (Ar)"^/^ 


+ 0 


+ 0 


(Ar)®/3 


1 

{A^ 

1 


vr^A^ 


1 - 


C2 


+ 0 


4(Ar)^/3 (Ar)^/^ 

and to the pressure ratio and anisotropy: 

a = 1 _ + o 

p± (Ar)2/3 ^ (At)4/3 ^ 


(Ar)^ 


(Ar) 


Ap 

e 


P± - P\\ 


= 2 


Cl 


+ Yya + o 


(Ar)2/3 (Ar)^/3 


(At) 


(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


6ln 


14 


the case 0 < s < 4 is studied. 


9 











































A is a parameter to be determined for each one of 


nn 1 ^ l + 21 og 2 

Stt 

the considered models. 

Understanding the results of the calculation of the solutions of the 5D Einstein equa¬ 
tions, hence, requires the computation of the boundary stress-energy tensor Ty in Eq. Q 
and the comparison with the proper time dependence following the previous equations. 
The various components of the stress-energy tensor of the boundary theory can be deter¬ 
mined using the holographic renormalization procedure developed in S ^ a nd shortly 
described in appendix [A} The functions ai{T) and 64 (t) in the expansions (18) and (20), 


computed through the Einstein equations, are required, and the energy density and the 
longitudinal and parallel pressures are obtained in terms of such functions: 


e — —^04 + €7 , 


p_L = -|a4 -h , 


P|| = 264 


(36) 

(37) 

(38) 


The functions e^yir), p_L, 7 (r) and P\\^^{t) are related to the distortion prohle 7 (r) in the 
boundary metric, and are reported in appendix [X| 

We can now discuss the three considered models, focusing on the effects of the bound¬ 
ary distortion and on the transient from the far-from-equilibrium state to the hydrody¬ 
namic regime. 


Model A. 

In the case of two distant pulses, the horizon is formed as a two-step process which 
follows the sequence of the boundary distortion, as one can infer considering the com¬ 
puted function A(r, r)/r^ which is depicted in Eig.[^ The outgoing radial null geodesics 
are clearly separated by a critical geodesic, the event horizon, which starts plunging back 
in the bulk immediately after each pulse in the boundary. This can be even more clearly 
seen in the plot of the r dependence of the effective temperature Tejj(r) and in the plot 
of the area of the apparent horizon per unit of rapidity, S(r/i,r)^, collected in Eig. 
The temperature has two decreasing regimes, the first one after the first pulse, with the 
system starting relaxation immediately after the end of the boundary distortion. Time 
intervals of about twice the duration of the pulse are very long with respect to the re¬ 
laxation time. The same phenomenon is observed considering the function S(r/i,r)^, 
which shows two time ranges in which it reaches a constant value, soon after the pulses. 
However, on such time scales, the components of the stress energy tensor have their own 
(computed) behavior, and no isotropization is observed after the first pulse. 

On the other hand, for a distortion represented by two close (nearly overlapping) 
pulses, no particular structure is observed before the end of the boundary perturbation. 
The horizon area per unit of rapidity, S(r/i, r)^, has a sharp increase during the distortion, 
and reaches the asymptotic large r value S(r/i,r)^ = tt^A^ soon after the end of the 
distortion. 


The hydrodynamic r—dependence of Te//(r) and of the energy density e(T), Eqs. (30) 
and (31), is reached immediately after the end of the distortion of the boundary r^, with 
values of the parameter A which are different for the two cases of different time intervals 
between the pulses. The values of A obtained from Te//(r) are collected in Table 
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Air,T]/r" 



Air,T]/r^ 



0 


-10 


-15 



Figure 1; Model A. The left panel shows the function A{r,T)fr"^, vs r and 1/r, obtained solving 
the Einstein equations Q-(IO) in the case of two distant (top panel) or close pulses (bottom panel) 
in the boundary metric. The color bar indicates the values of the function. The gray lines are radial 
null outgoing geodesics, the dashed dark blue line is the apparent horizon, and the continuous cyan 
line is the event horizon obtained as the critical geodesic. The excision in the low-r region, used 
in the calculation, is shown. In the right panels, the 3d representation of j4(r, r)/r^ is displayed 
together with the apparent horizon (continuous line). 


The temperature can be obtained from the condition S(r/i,r) = 0, or from the energy 


density Eq.(29), with numerical results differing at the level of 2%. For a quantitative 
discussion, a time t* can be fixed imposing that the energy density e(r*) differs from 


the hydrodynamic value eniT*) in (31) by less than 5%: 


e(r*) - enir*) 


€{t*) 


= 0.05. 


(39) 


In both cases of distant or overlapping pulses, t* essentially coincides with the end of 
the boundary metric distortion. From the plot of Ty in Fig. we observe that at 
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T* the pressure anisotropy is still sizeable, and the values of Ap/e and p\\/pi_ are still 
much larger than the ones expected by viscous hydrodynamics. Therefore, we define the 
’’thermalization time” Tp from the condition 


P\\{Tp)/Pl{Tp) - {p\\{Tp)/p^{Tp)) 


H 


P\\{tp)/pl{tp) 


= 0.05 


(40) 


where {p\\{Tp) /p_\_{tp))h is given by (34). Quantitatively, from the condition that p\\/p±_ 
differs by less than 5% from the asymptotic NNLO expression, we can set Tp = 6.8 for 
two distant pulses. The difference Tp — t* can be expressed in physical units if one scale 
in the system is set. Imposing T(,ff(T*) = 500 MeV corresponds to Tp — t* 0.60 fm/c, 
which indicates the elapsed time between the end of the pulse and the restoration of 
the hydrodynamical regime. For two overlapping pulses, the time difference is — r* ~ 
1.03 fm/c. The numerical values are collected in Table Notice that r* and Tp differ in 
general, although Eqs. (27) and (28) hold. Indeed, at a generic r the relation between 
energy density and its asymptotic limit can be written as: 


e(r) = ^(r)eH(r) 


(41) 


with ^(r*) such that (39) is verified. For, e.g., the longitudinal pressure, using (27) and 
(41) one has P\\{t*) = C{t*) P\\,h{t*) — Therefore, due to the last term, 

the value Tp verifying (40) does not necessarily coincide with r*. 


Model B. 

The purpose of this model is to investigate a boundary distortion characterized by 
two largely different time scales. The horizon is formed immediately after the boundary 
deformation is switched on, as depicted in Fig. then starts plugging down in the bulk 
until the short pulse becomes active. This can be seen clearly in Fig.[^ where the effective 
temperature and display a decreasing and constant regime, respectively, in the proper 
time interval r* < r < roq, after which a violent perturbation takes the system out of 
equilibrium. The time t* can be fixed using (31). For t > t* Tt t^ the determination 
of A from ^e// and the energy density e gives A = 1.12, with a variation of 3% around 
this value. Thermalization is reached at the time Tp obtained by the condition (40), as 
one can infer considering the observables depicted in Figs. and Setting the scale so 
that the temperature is Teff{T*) = 500 MeV, we hnd Tp — t* 0.42 fm/c. 


Model C. 

This model is constructed with the aim of representing a situation close to the physical 
system during the scattering process, with two different time scales and several pulses of 
various number and intensity. Immediately after switching on the boundary deformation, 
the horizon is formed in the bulk geometry and follows the profile of the distortion. Even 
in short time intervals between the pulses, the horizon starts to plung in the bulk, with 
a decreasing temperature and a saturation of This can be inferred considering the 
computed function A/r^ depicted in Fig. and, in details, studying the r dependence of 
the effective temperature and horizon area per unit rapidity in Fig.[^ S^(r/i,r) displays 
a step-like behavior in proper time, closely following the distortion 7 (r), and reaches a 
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Figure 2: Model A. The panels, for the prohle 7 (t) with different time intervals between the 
pulses (left and right columns), show (from top to bottom) the function 7 (r), temperature Te//(r), 
horizon area S^(r/i,r) per unit of rapidity, and the three components e(r), p±{t) and p\\{t) of the 

stress-energy tensor Ty (apart from the factor ). 


constant value from t = Tj on. On the contrary, the various components of the stress- 
energy tensor have structures which become regular only for r > rj?, after the last pulse. 
As one can observe analyzing the results in Fig. the energy density reaches the NNLO 
hydrodynamic behavior at r = r* — Tf, while pressure anisotropy persists for longer 
times. The parameter A, obtained from different observables, takes the value A = 1.59, 
with a variation similar to the one in models A and B. The pressure anisotropy Ap/e 
and ratio p\\/p± set the value Tp of proper time, as one can infer considering the results 
in Fig. setting Teff{T*) = 500 MeV, we obtain Tp — t* 0.2 fm/c. 
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Figure 3: Model A. From top to bottom: temperature Teff{T), components e(r), p±{t) and 
of the stress-energy tensor, pressure anisotropy Ap/e = {p± — p\\)/e and ratio p\\/p±, computed for 
T > for the boundary distortion with two distant (left) and two overlapping pulses (right panels). 
The short and long dashed lines correspond to the hydrodynamic result and to the NNLO result in 
the 1/r expansion. 

5 Conclusions 

We have investigated the effects of different types of distortions of the boundary metric, 
in a boost-invariant setup. The quenches are introduced as a way to take the boundary 
theory out-of-equilibrium, and the distortion profiles are used to describe processes with 
different time scales and intensities. As a common feature, we observe a rapid formation 
of the horizon in the bulk metric, with the possibility of defining an effective r—dependent 
temperature. In coincidence with the end of each main distortion, the relaxation starts 
with the horizon plunging in the bulk. We have seen how sequences of distortions break 
the relaxation process. We find that for all the considered distortion profiles, Teff{T) 
starts to follow a viscous hydrodynamic expression as soon as the quench is switched 
off. At this time the pressures are different, and evolve towards a common value which 
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A[r,T]/i^ 



Figure 4: Model B. Function A{r, T)jr‘^, vs r and 1/r, outgoing radial null geodesics, event horizon 
and apparent horizon. Lines and colours have the same meaning as in Fig. 


model 

T* 

Tp 

A 

At = Tp — t* (fm/c) 

A(l) 

5.25 

6.8 

2.25 

0.60 

A (2) 

3.25 

6.0 

1.73 

1.03 

B 

5 

6.74 

1.12 

0.42 

C 

9.45 

10.24 

1.59 

0.20 


Table 1: Numerical values of the relevant quantities for models A, B and C of boundary 
distortion. 

is reached at a later time. Setting = 500 MeV, the elapsed time before the 

restoration of the hydrodynamic regime is always a fraction of fm/c. Although the 
system described by the holographic approach is in several respects different from the 
real QCD system, and boundary sourcing a quite abstract representation of the heavy 
ion collision process, the obtained results allow us to argue what can be expected in 
realistic situations. 
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A Boundary stress-energy tensor 


For the sake of completeness, we briefly describe the main steps of the calculation of the 


boundary stress-energy tensor leading to Eqs.(36)-(38). The bD metric (14) can be 
written in ADM form (3^, 


ds^ = N^dA + {dy^ + M^dr) {dy^ + M'^dr) 


(42) 
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Figure 5: Model B. The panels show (from top to bottom) the profile 7 (t), temperature Teff{T), 
horizon area per unit of rapidity Y?{rh-,T), and the three components e(r), p±{t) and of TH'. 


in terms of the indnced metric of the lapse function N and of the shift vector field 
The stress tensor for the space foliation M.r obtained at constant r, is given by 




6S 




flU 


(43) 


with h = det [h^y) and S the gravitational action. The boundary stress energy tensor is 
obtained from the limit: 

(44) 


= lim 


The action S in (43) includes as a counterterm a local functional Set of whose 

. As a consequence, 


contribution to regularizes the divergences at r 
the stress-energy tensor is written as 


oo 


33 


31 


/yS / 1 1 \ 

= lim ^ - 3/1^^^ + - logr , 

r^oo y 2 J 


(45) 
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Figure 6: Model B. From top to bottom: temperature Teff{T), components e(r), p±{t) and 
of the stress-energy tensor, pressure anisotropy Ap/e = {p± — p\\)/e and ratio p\\/p±, computed for 
r > Tj. The dashed lines correspond to the NNLO expressions. 


with the extrinsic curvature, the boundary Einstein tensor with curvature 

tensors defined with respect to h^y. The counterterms ^ cancel the 


-2 


powers from to r~^ from the first two terms in (45), and introduce terms of order r 
which contribute to the final result (45) 33 . The last term proportional to , with in 
our case 


afj,y = diagl --a 4 (r),e' 


o7(r) 


-^«4(t) +2/?4(r) 


=7('r) 


-2“4(r) + 2/34(r) 


,e-27(r)^2 


--a4(r) - 4/34(r) 


(46) 


and a 4 (r) and j5^{T) in (15)-(17), cancels the ^logr contributions in (45). The results 
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A[r,r]/r2 



Figure 7: Model C. Function A{r, T^jr^ ^ vs r and 1/r, outgoing radial null geodesics, events horizon 
and apparent horizon. The lines and colors are indicated as in Fig. 


for the components of the stress-energy tensor Eq. 0 are in Eqs.(p^-(p8l), with: 


6^ - 


25 y , ^ // 19 /3 , Ty/.^iy/ _ y^(3) 

^07 96 / w 144 / , .-jO / ^ 96 ' V M 48 / 


144 7-3 


+ 


_ _ 1 //2 , Av ^( 3 ) 

256 ' 8 ' 32 ' ' 


(47) 


1 40 q V , 145 //N// _ J_ (3) 

~ __ t _I 7 _ 9 > ^ 1728 ^ 4_ 192 ' ^ 288 ' v / / 288 ' 

6r4^ 576r3 r2 ^ r 

_ ily4 _ ^y2y' _ J_ "2 , Ay^(3) , J_^A) 

256 ' 64 ' ' 96 ' 96 ' ' 192 ' ’ 


(48) 


P\\ 7 


QKO y 49y2 437 // Ay3 2i« / // x^(3) 

-Joy 7 , 36 / 4^ 864 ' 12 ' ^ 288 ' ' ^ 72 ' 


3r4 288 r3 


_ ily4 . ^y2y' _ ^y'2 , Ay^(3) _ 

256 ^ 32 ' ' 96 ' 96 ' ' 96 ' 


(49) 


Energy density and pressures are hence affected by the variations of 7 (r). The above 
expressions coincide (after the correction of an overall minus sign in e^) with the ones 
in 12 after the recognition of terms due to the regularization scheme proportional to 
a^u, related to the coefficient of the log[r] term in the large-r expansion of the metric. 


The covariant conservation of the stress-energy tensor gives again Eq. (21) 
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Figure 8: Model C. From top to bottom: profile 7 (t), temperature Teff{T), horizon area per unit 
of rapidity S^(r/i,r), and the three components e(r), p±{t) and of Ty . 


B Testing the numerical algorithm 


To quantify the accuracy of our numerical algorithm, based on the Runge-Kutta method 
for the solution of the differential equations in the variable r, we monitor the ratios 


Ti%{r, t) 


7^7(r, r) 


T^sir, t) 


S(S)' + 2S'S-2S2 
|S(S)'| +2|S'S| +2|S2| 

S(R)' + I (r.'B + B't'j 
|S(Ry| + §(|S'R| + |R'S|) 

A' + ‘iB'B - 12^ + 4 
\A''\ +2,\B'B\ + |12^| +4 


(50) 

(51) 

(52) 
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Figure 9: Model C. From top to bottom: temperature Teff{T), components e(r), p±{t) and 
of the stress-energy tensor, pressure anisotropy Ap/e = {p± — P||)/e and ratio p\\/p±, computed for 
r > T^. The curves have the same meaning as in Fig. 


7 ^9(r, r) 


S + i 


1 

\n+i{ 

J52S| + 


T^ioir,T) = 


S" + 

|S"| + |S'2|S| 


(53) 

(54) 


in the (r, r) domain behind the excision r > rmin- The final results are shown in Fig. 10 
in the case of model 13. The resulting ratios TZq, TZj, TZg deviate from zero at a level 
smaller than 0(10““^), but for a tiny region close to the excision in the case of TZq and 
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Figure 10: Ratios 7^j(r, r) in the case of model B. The white region corresponds to the excision in 
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0.04 




Figure 11: Ratio Tlg{r,T - 
curve) and dr = 0.5 • 10“^ 



= 4) computed in the range 1.2 ^ r ^ 3.0 for dr = 10 ^ 
(purple, dashed curve) in the case of model B . 


(pink, continuous 


TZs, and for a few spots in the range 3 ^ r ^ 5 (where the source function 7 (r) is peaked) 
in the case of TZ^. A similar result is obtained for TZg and T^io, with larger deviations 
from zero: this is a consequence, for TZiq, of the smallness of the two addendi in (10), and 
for TZg of the computation of S = clrS + \A (S)' through a discretized time derivative. 
To monitor the time slicing, we compute TZg for two different time steps in the region 
1.2 ^ r ^ 3.0 at r = 4, where there is the largest deviation from zero. Comparing 
the results for dr = 10“^ and dr = 0.5 • 10“^, as shown in Fig. 11 one observes the 
decreasing of TZg by doubling the grid points. 
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